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Abstract 

O 

^p -1 , We present general exact solutions for two classes of exponential potentials 

in a scalar field model for quintessence. The coupling is minimal and we 
' consider only dust and scalar field. To some extent, it is possible to reproduce 

experimental results from supernovae. 
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I. INTRODUCTION 

In the last few years, new models of the universe have been built taking dark energy 
into account JT|^] . Together with baryons, cold dark matter, photons and neutrinos, a fifth 
component has been added, the socalled quintessence field Q ( or ; i n general, the x-field 
Generalizing ideas like that of a cosmic equation of state variable with a A-term | |10| . 



with respect to a more usual cosmological constant A, such a Q-field, even if still implying 
a negative pressure contribution to the total pressure of the cosmic fluid, is characterized 
by the fact that its equation of state is given by —1 < wq = Pq/pq < 0, pq and pg 
being, respectively, the pressure and energy density of the Q-field. Actually, the interval 



— 1 < wq < —0.6 is usually considered |TT|. As a matter of fact, when wq = — 1 we recover 
a constant A-term |T2|-|14|, which can be regarded as a measure of vacuum energy density, 
leading to the well known discrepancy between theory and observations [|T3|,|T2|,|T^| , based on 



the question of why pq is so small with respect to typical particle physics scales. (But there 
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are also mechanisms of relaxation of the cosmological constant during the initial inflationary 
stage, which could explain such a discrepancy; see [fL6|| , for instance.) 

As well known, an interesting possibility to handle the presence of quintessence in the 
universe is to see it as given by a scalar field ip slowly rolling down its potential V(<p). If we 
define 

P V = \^ + V{^) ■ p^ = ^ 2 -V( V ), (1) 

(dot indicating time derivative, and V(ip) being the potential for ip), the slow rolling condi- 
tion immediately gives p v < and w v = p v / p v ~ —1. With such a negative pressure, the 
universe evolves like in a sort of present day soft inflationary scenario, so allowing to explain 
observations on supernovae |17|-fjl]] and why vacuum and matter densities are today compa- 
rable ('cosmic coincidence' problem |2^JT^Ji3|] ) . (Alternatively to quintessence, a negative 
pressure and an explanation of current observations can also be obtained in a Chaplygin 
cosmology [p3fl .) 

Many cosmological models with a dynamical scalar field have been proposed, showing 
scaling solutions, i.e., such that at some time p m and p v simultaneously depend on some 
powers ni and n 2 of the scale factor a, acting as attractors in the phase space. When ni = n 2) 



we have the socalled self-tuning solutions [24|, which are typically driven by exponential 



potentials. This kind of potential has been studied extensively |25|-p4|, especially from a 
qualitative point of view (see [|35| and references therein, for instance). 



The simplest possibility is of course V((p) = ae^, which is often discarded (see discussion 
in |6|,[14]]). More promising seems to be a combination of two terms V(<p) = ae X(p + (3e~ Xip 



32] 



In this paper, we consider a particular class of both these types from a different point 
of view, obtaining general exact solutions. This allows a very stringent comparison with 
experimental data on supernovae, so that also the first type seems to deserve further investi- 
gations; for the second type, we obtain a solution which can mimic very well the presence of 
a cosmological constant in the late evolution of the universe. Both of them are not, strictly 
speaking, scaling solutions, although this concept may be recovered in a more general sense. 

Another experimental fact which we use is the strong evidence of a spatially flat universe 
||36|| . Thus, we set the scalar curvature k = in all our equations from the very beginning. 
However, as we shall see, the values of Q m and Q\ (Q(p in our case) derived from the 
experiments strongly depend on the model, so that some discussion is needed. 

Mostly, the scalar field <p has been considered as minimally coupled to gravity, even 
if (more recently) a nonminimal coupling has also been introduced P71-S3]. Here, we will 
consider a very simple model consisting of a two-component cosmological fluid: matter and 
scalar field. 'Matter' means baryonic + cold dark matter, with no pressure, and the scalar 
field is minimally coupled and noninteracting with matter. Clearly, this model cannot be 
used from the very beginning of the universe, but only since decoupling of radiation and dust. 
Thus, we do not take into account inflation, creation of matter, nucleosynthesis, etc. The 
main shortcut is that we cannot really check for the tracker feature |5]|| of the ^-solution. 

In Sec. II we take a particular exponential potential into account, deriving the general 
exact solution of the cosmological equations and, thereby, cosmological parameters as func- 
tions of time, so allowing the comparison with observational data. Sec. Ill is devoted to 
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the same kind of considerations for a potential given by a linear combination of two such 
exponentials. Conclusions are drawn in Sec. IV. 



II. AN EXPONENTIAL POTENTIAL 

A. Mathematical treatment 

Let us consider a spatially flat, homogeneous and isotropic universe, filled with two 
noninteracting components only, i.e., pressureless matter (or dust) and a scalar field ip, 
minimally coupled to gravity. The cosmological equations are then 

3H 2 = ^f(p m + p v ), (2) 
H + H 2 = -^(p m + p tp + 3( Pm +p (p )) j (3) 

(p + mcp + v'{ip) = o, (4) 

where prime indicates derivative with respect to (p, H = a/a is the Hubble parameter, 
p m = w m p m and p^ = w^p^p are the equations of state for matter and scalar field. Let us 
stress that is not constant, and that we want to describe some features of cosmology 
after the decoupling. 

We set w m = 0, so that p m = Da~ 3 . The parameter D = p m0 a 3 (the lower index '0' 
indicating present day values) is the amount of matter. The equations can also be rewritten 
as 

+ 3 + V'foO = 0. (7) 

In this Section, we consider the potential 

V( V ) = B 2 e-^, (8) 
where B 2 is a generic positive constant and 

a 2 = (9) 

(The minus sign in the exponential is irrelevant, since there is symmetry with respect to a 
change ip — > —(p.) 

This type of potential leads to a late time attractor in a scalar-field dominated situation 
(Q v = l,w<p — —0.5) [|rj],|32|. Being aware of such a behaviour, anyway, we stress that we 
are especially interested in the contemporary or, at most, the recent past regimes, where 
the situation is different. Usually, associated with an exponential potential, a scalar field 
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is considered such that ft v = 8nGp ip /(3c 2 H 2 ) is practically constant during part of the 
matter-dominated era. This implies that assuming w v ~ constant leads to a constant 
ratio of quintessence to matter energy density so that Q v (being < 0.15 at the beginning 
of matter-dominated era, due to nucleosynthesis |]29|j30[| ) must remain small forever |6"1,|T4 



Mainly for such reasons, this kind of potential is not considered as suitable for a quintessence 
field. 

The particular choice of Eq. (|9|) for o allows for general exact integration of equations. 
Such a choice was in fact used in the context of inflationary theory by us and others 

5B], with a scalar field only. 



Let us concentrate on the second order equations (6) and (7), while Eq. (5), which is 
a first integral, is considered as a constraint on the integration constants. Let us introduce 
the new variables u and v, defined by the tranformation 

1 u 

a 3 = uv ; <f — log—, (10) 

a v 

which is always invertible (the Jacobian being J = 2/cr). We get for the potential 

V(u,v) = B 2 ^, (11) 

and Eqs. (6) and ( 7) become 

u = ; v = luu, (12) 
where uj = o~ 2 B 2 = 12irGB 2 jc 2 > 0. They are immediately integrated to 

u(t) = uit + u 2 , (13) 

v(t) = \umjt 3 + \u 2 ut 2 + vxt + v 2 , (14) 
o 2 

being Ui, u 2 , t>i, and v 2 arbitrary integration constants. Taking into account Eq. (5), we 
find 

2 1 47rG ^ n 

-u x v x - -uju 2 —D = 0. (15) 

o o C 

Since D is a physical parameter, it would be natural to use it as given and derive one 
of the other constants. But this complicates calculations without substantial advantages, as 
its value depends on the normalization of the scale factor a, which can be fixed arbitrarily. 
Thus, we will determine D from Eq. (|15|) 

c 2 

D = (2u x vi - uju 2 ) . (16) 

lZ7rCx 

Being of course D > 0, this gives limitations on the choice of u\, u 2 , V\. 



The well known [^,0] solution <p — 2/ a/3(1 + a) logt, coming from the potential 
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is a very particular case of what we find in Eqs. (|T3|) and (|14| ). It can be obtained by setting 

B2 = on7 a l ; °= \/3(l + «) ; u 2 = « 2 = v 1 = 0. (18) 
o(l + oj) z 

Eq. (PS|) then gives D = 0, so that we get a model without matter, not really interesting 
in our context. 

A more interesting possibility is given by the choice u\ — V\ — v 2 = 0, involving, from 
Eq. (|I~6D , u 2 < 0. This implies, in fact, 

„ c 2 ujuo □ 1 9 9 4 . , 

= - 127rG' ^=2"^. = < 19 » 

In this particular case, (oc a -3 ) scales as p m . (But see below for a discussion on the 
scaling properties.) 

Without any special assumptions on constants, we can get many important quantities 
as functions of u and v (we do not write them explicitly in terms of t, for sake of brevity) 

1 / (uv — uv) 2 u\ , , 

P * M = ^{ L ^ L+ "v)' (20) 
, . 1 / (uv — uv) 2 u\ ,^ . 

*fc" ) = ^(4=w u s)- (21) 

(Mt> — Mt>J z + ZtUU^f 
rT / x UV + UV ,„„ s 

=^ , 23 

n v (u,u) = ^ ' , 24 

r i 111) J- -7/-7J 1^ 



24ttGDuv Auuvi) — 2uu 3 v 
ttm{u, v) = = — — — ■ — — , (25) 

C Z [UV + UV) Z (UV + uv) 



the last equality coming from Eq. flIB|). It can be easily checked that f2 m + ft v = 1. 



B. Physical considerations 

We now pass to simplify the situation with a suitable choice of initial conditions. Of 
course, the following choice is not the only one possible. We shall see, anyhow, that it is 
able to reproduce observational data. Other more general choices may improve the situation, 
but we do not treat them here. 

If ti n = —u 2 /ui (which is always possible, being u\ ^ 0), the scale factor a is zero, and 
we can show that there is no other time t > t in when this occurs again. We can thus fix 
the time origin in such a way that a(0) = 0. This condition has to be interpreted just as 
an arbitrary choice of the time origin. The real beginning (of physical meaning) for the 
model starts a little bit afterwards, at a time t\. This delay is otherwise arbitrary, so that 
this setting does not seem to exclude important CclSGS, clS said before, and leads to a great 
simplification in the formulae. Now, a(0) = implies u 2 = or v 2 = 0, or both. If we 
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set only one of them to zero, we obtain <p(0) = oo (which could be accepted, but is rather 
disturbing), and, most of all, Q v (0) = 1, which would mean an initial scalar-field dominated 
universe, with a neglectable content of other types of matter. Of course, if we consider 
the situation in general, the scalar field does dominate. But, as already mentioned, in our 
case we start after decoupling time, when a matter-dominated behaviour seems to be more 
natural. So, if we set U2 = v 2 = 0, we get instead (in a matter-dominated situation, then) 

I u 1 c 2 
<p(Q) = — log^ ; ft„(0) = ; D = — u 1 v 1 . (26) 
a V\ 07rG 

We prefer to stick to this choice, so that we have 

u(t) = uit ; v(t) = ^uiut 3 + v 1 t. (27) 

6 

Let us now define a time scale t s such that H(t s ) = l/t s , which is of the order of the age 
of the universe. This leads to consider a dimensionless time r = t/t s . From Eqs. ( p3|) and 
(0) we get 

t = — . 28 

UJUi 

By means of these choices the formulae found above for the relevant cosmological pa- 
rameters reduce to 

Pf = tr 2 t ( 2( 1 + T 2)2' ( 29 ) 

2(3 + 2r 2 ) 

Pip ~ " ff 2 t 2 (l + r 2 ) 2 ' (30) 
3 + 2r 2 

W - = ~3T4^' (31) 

a = (u l v l t 2 s (l + T 2 )T 2 f\ (32) 

(1 + 2)3 = ^TT7 2 T' (33) 

H ~ 3t s r(l + r 2 )' (34) 
1 +r 2 

fim = (l + 2r 2 ) 2 ' (35) 

where z = a(ro)/a(r) — 1 is the redshift, and tq indicates the present time. 
If we define dimensionless pressure and energy density 

0H 2 0H 2 

p<p = -^-Pf ; = ~y~p^ ( 36 ) 

we find the equation of state for the scalar field 



p v = p v -\2 + Gy/A-p v , (37) 



6 



which is well approximated by 



^ = -0.382^-0.196^;. (38) 

In Fig. 1, we compare the plots of the two functions in Eqs. ( E571) and ([53), and show 
that the approximation is quite good; a comparison is also made with a straight line = 
—0.86/5^,, where the coefficient has been obtained through a numerical approximation, as 
well those in Eq. (|3~8D . 

From Eqs. ( |2"9"D and (^) it is possible to derive 

_ dlogp^ 3r 2 (5 + r 2 ) 

l ~ dhga (3 + r 2 )(l + 2r 2 )' 1 ' 



Thus, it is clearly 



Si — > for r — ► 0, Si — > -- for r — > oo, (40) 



so that we asymptotically have two scaling regimes: p v ~ const, for early times, and oc 
a -3/2 f or | a ^- e times. j n f ac t ; this approximation holds for a very long time, well far from the 
asymptotic values. Indeed, computing the n-th derivative S n = a(d5 n ^i/ d(r 2 )) / (da/ d{r 2 )), 
it is possible to show that they are all asymptotically zero up to, say, n = 10. 

An estimate of the moment in which the regime changes can be given finding the max- 
imum (or the minimum) of 62- This is achieved at r ~ 0.3, and it is remarkable that this 
result depends only on t s . All this discussion shows that it is possible to generalise the 
concept of scaling solutions. The situation is illustrated in Figs. (2) and (3). 

As a matter of fact, in the literature it is widely accepted that using an exponential 
potential leads to a dark energy density which scales like matter. Our results seem to be 
in contrast with this statement, which is a consequence of assuming almost perfectly 
constant. It is not our case, as shown in the following. 

As in ||, we can use the function T = V"V / (V) 2 . Defining 



(fi 2 l+w v . dhgx .. d 2 \ogx . 

X = ^T7 = 1 ' x = 77i ' x = 77i 2> 41 

2V 1 — w tp d log a d log a z 



it is possible to find 



r _ 1 w m -w ip l + w m - 2w v x 2 x 

2(1 + w v ) 2(1 + w v ) Q + x l + w^(Q + x) 2 ' [ ' 

If one makes the assumption that w v is nearly constant, then x ~ x, ~ 0. Since it is T — 1 
strictly in the case of our potential, Eq. (|4"2"| ) implies w m ~ w v . The point is that it is not 
true that our exact solutions for the exponential potential lead to x ~ x, ~ 0. For instance, 
we have 

r 2 3 

x = —, — , x = , (43) 

3(1 + r 2 )' l + 2r 2 ' V ' 

so that x — > only asymptotically. Since in our model r < 1, x and x are then far from 
being zero. 
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But let us also consider the third term in the right-hand side of Eq. ( f4^ ) with w r , 
(matter is simply dust in our model). Substituting our solution, we get 

1 + w m — 2w v x 9 + 8r 2 



2(1 + w v ) 6 + x 12r 2 + 16r 4 



(44) 



we see that this expression diverges for r — > 0, being always > 0.5 in the useful range of r. 

In our opinion, the main check for the solution in Eq. (|27|) is thus only its capability to 
reproduce the experimental results, which we are going to do just below. 

From Eq. (|35|) we get 



To = ^ • (45) 

Once we give an acceptable value for O m0 , we obtain a value for To. For instance, Q m0 = 
0.3 gives r = 0.82, and this implies w vo = —0.76. If the value H = 100/i km s^Mpc' 1 is 
also given, we get 

2(1 + 27?) = L14 
* s 3^(1 + r 2 ) flo' W 

assuming /i = 0.7, we have t s = 15.8xl0 9 years, and to — 13xl0 9 years. It is also possible 
to obtain the relation between Wtp and Q m0 



w 



1 + 8fi m0 — 3-y/l + 8fi m0 
vo ~ 4(1 - fi m0 ) • 



(47) 



For Q m0 = 0.2 4- 0.4, we get w vo = —0. 699 H — 0. 811, and the value —0.5 is reached only in 
the case of f2 m0 = 0- It is also possible to obtain function of the redshift 

w <p = , = == , (48) 

where £ = (1 + z) 3 - 

Figs. 4, 5, and 6 show the dependance of w — (ft versus r, fi m0 , and (, respectively. We 
can see that varies much and has values ~ —0.78. For ( = 1 (now) it is w v = —0.77, 
and already for ( = 4 (z ~ 0.59) we find w v ~ —0.89. It is remarkable that the values of 
\w<p\ are greater than 0.7. According to our knowledge, this feature is found only in p4 |. 



Another interesting quantity is the present value of <p. After staightforward algebra we 

get 



a 



l l ° g " 4a "° + A/l + Sfimo)) > (49) 



and we see that this value depends on the observed parameters and on the value of B 2 , 
which was until now completely undetermined. Now, for r = 0.82 and to = 13xl0 9 years, 
we have B 2 = 2.5x 10 -47 exp(— ery? ) GeV 4 . Considering ip ~ 1/6 M P (M P being the Planck 
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mass) we see that exp(— atpo) ~ 1, and we can determine the unknown parameter for the 
potential 

V{0) = B 2 ^ 2.5x 10~ 47 GeV\ (50) 

But we have also to observe that a 'little' change in ip entails a 'large' change in B 2 . 
For instance, if <^ ~ Mp, then exp(— o"<^ ) ~ 0.0025 and B 2 changes of three orders of 
magnitude. Due to Eq. (|26|), we have that ip(0) = —l/<T\og(2irGB 2 t 2 /c 2 ). This means 
therefore that a relatively wide range of initial values of ip ends up to a narrower set of final 

Thus, everything seems to work fine, but things are more complicated. Indeed, one has 
to ask what is really measured in the supernovae experiment. The value f2 m0 = 0.3 is not 
a direct consequence of the data, since it depends on the model, which uses the constant 
A-term. What we really measure is the distance modulus, so that it is this quantity that 
we should compare in the two situations. Here, we limit ourselves to a very qualitative 
discussion. 

Let us recall, then, the definitions of luminosity distance (in Mpc) 

^ = 3000(1 + ^) J — , (51) 

and distance modulus 



5 = m — M = 5 log 10 cIl(z) + 25. (52) 
We have thus to compare this last quantity in the case when H(z) is taken from the 



usual model with A fl8|j21 |, that is 



H{z) = H oy /(l + ^) 2 (1 + <W) - z(2 + z)(l- n m0 ), (53) 

with the one obtained eliminating r from Eqs. (|33|) and (|34|). 

In Fig. 7 we compare 5 with 5 (let us mark with a ~ the values for the model with 
A). The agreement is almost perfect, up to 0.06%. But there is a trick! 5 was obtained 
from a value Q m o = 0.37 . Of course, this value is still in the possible range but at its 
limit. If we decide to trust strongly on the value Q m Q = 0.3 and want to obtain the same 
good agreement, we have to change the value of r to 1.22. This gives a very different value 
fl m0 = 0.16 in the model with (p. This is again at the limit of possible estimates (due to 
other investigations on dark matter). 

In conclusion, we see that this solution (with the potential in Eq. (||)) is indeed difficult 
to fully adapt to observed data, but for reasons which are not easy to investigate without 
general exact solutions. Moreover, it is not clearly incompatible (until we get better data); 
therefore, it seemed to us useful to present it in detail. 



III. TWO EXPONENTIALS COMBINED 
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A. Mathematical treatment 



We now consider a combination of two exponentials, which will give us much better 
results, as expected. The procedure strictly follows the above one. 
Let us consider the potential 

V(ip) = A 2 e a ^ + B 2 e~ aip , (54) 

with a 2 = 12nG/c 2 as before, and A 2 , B 2 arbitrary parameters. We use, now, the following 
change of variables 

3 u 2 -v 2 1 B(u + v) 

a = — - A — ; y? = -iog^7 r, 55 

4 a A(u — v) 

which is invertible, provided that a^O. This leads to 

u 2 -\- v 2 

V(u,v) = 2AB- -. (56) 

With these variables Eqs. (6) and (7) are rewritten as 

u = oj 2 u ; v = —u 2 v, (57) 

where now 

2 12nGAB 

uj 2 = - 2 . (58) 

Again, the integration is immediate, and gives the general solutions 

u(t) = ae ut + Pe~"\ (59) 
v(t) = v 1 sin(cjt + v 2 )i (60) 

with a, (3, vi, v 2 arbitrary constants. As before, we derive D from the constraint in Eq. (5) 

c 2 u 2 {v l + Aaf3) 

24ttG ' 1 ' 

Being D > 0, this implies v\ < —4aj3. A change in the sign of v\ has the only effect of 
changing the sign of tp, and interchanging A 2 with B 2 . So, we can set V\ > without any 
loss of generality (the case v± — is obviously equivalent to considering a A-term). As a 
consequence, a and (3 must be non zero and with opposite signs. 

Again, we can write some important functions in terms of u, v 

2{{vu-uv) 2 + uj 2 {u'-v*)) 
Pf~ a 2( u 2_ v 2y ' W 

2((vu - iivf - u 2 (u 4 - v 4 )) 

(vu — uv ) 2 — u 2 (u 4 — v 4 ) 
w v = 7 ro 57—3 tt- 64 
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(This last expression gives w v ~ —1 when (vu — iiv ) 2 < a 4 - t> 4 , which certainly can 
happen for sufficiently large times.) 
For the redshift we have 



and finally 



An 3 

:i+-) 3 = ^4, (65) 

u 2 — V 2 



— — i . 2{uu — vv) ,„„ s 

H M = in2 — 2V' 66 

3(u 2 — v 2 ) 

(u 2 -v 2 )(v 2 -u 2 + uj 2 (u 2 + v 2 )) 

U m (u,v) = — — , (67) 

(WH — VV) 

^ / x (vu — iiv) 2 + uj 2 (u 4 — v A ) , . 

v) = y - V ) 6g 

(uu — vv) 2 



with Q m + fltp — 1, of course. 



B. Physical considerations 

The potential in Eq. ([54]) has a nonzero minimum V m i n = A 2 + B 2 . This is unusual in 
quintessence theory, since V m i n = is optimal to remove the fine-tuning problem. Indeed, 
Vmin 7^ can be seen as a disguised A-term. Anyway, in the following we find that this 
is actually the case when the scalar field is almost stationary near the minimum of the 
potential. But there is also the possibility of a slower rolling far from the minimum. The 
situation is then similar to that in Sec. II, but the additional term in the potential now 
allows to achieve a better agreement with observational data. 

Let us now make a trial for the choice of the free parameters. We set again a(0) = 0, and 
ask for nonsingular ip(0). The situation, and hence the interpretation, is the same as above. 
Thus, we pose a = —f3 = A/2, t> 2 = 0. It is also possible to fix an arbitrary normalization 
for a and set V\ — 1, obtaining at last 

u(t) = Asinh(o;t) ; v(t) = sm.(ut). (69) 

We get now 



implying |A| > 1, and 



We define a dimensionless time r = cot and get 
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. . /A 2 sinh 2 r - sin 2 r\ 1/3 ,„ n s 
a(r) = ^ J , (72) 

. o A 2 sinh 2 r — sin 2 t . , 

(l + ^) 3 = X2 ■ , 2 • 2 ' (73) 

\ z sinh t — sin r 

^(sin(2r)-A 2 sinh(2r)) 
3(sm r — A z sinh r j 
^ ^ A 2 (cosh r sin r — cos r sinh r) 2 + (sin 4 r — A 4 sinh 4 r) 
A 2 (cosh r sin r — cost sinh r) 2 — (sin 4 r — A 4 sinh 4 r)' 
2(A 2 - l)(coB(2r) + A 2 cosh(2r) - 1 - A 2 ) 
ml J (sin(2r)-A 2 sinh(2r)) 2 ' 1 } 

In comparison with the situation in Sec. II, we now have one more free parameter. This 
gives the possibility of a much better agreement with observational data. 

We can also take into some account the final and initial values of ip, i.e., <p and ipi = y?(0). 
If t is the present dimensionless time, we get 

/ / xx A 2 sinh To + A(sinr — sinhr ) — sinr 

exp (a(ifii - <po)) = , — — — : . (77) 

A z smhr — A (sin To — smn To) — sin To 

When A >> 1, we have ipi ~ v?o • In fact, tp is practically constant and we have — — 1, 
with nearly perfect emulation of a cosmological constant. On the other hand, if A ~ 1, then 
exp (a(ifii — ifo)) ~ 0. This can be interpreted as ifii ~ oo, or, better, as the possibility of a 
wide range of </?j's, with nearly the same final ip . 

Whatever is A, anyhow, it is possible to obtain a good agreement with observational 
data. We here give only two extreme cases, with A = 30 and A = 1.1. 

In the first case (case I) we set 

A = 30 ; t = 1.2 ; u = 2.8x lO^V 1 = 2.8x 10~ 42 GeV, (78) 

which gives 

#o = 70 ; fi m0 = 0.3 ; w vo = -0.999. (79) 
In the second case (case II) posing 

A = 1.1 ; T = 0.44 ; uj = 1.07xl0" 18 s _1 = 1.07x 10~ 42 Ge\/ (80) 

gives 

t7 = 70 ; n m0 = 0.3 ; w VQ = -0.76. (81) 

Again, we find that the value of w v is less than —0.7, as in Sec. II. In Figs. 8 and 9 
the distance modulus 5 for these two cases is compared with the A-term case. As before, 
the agreement is quite good, but now we have Q m0 = 0.30, so that, with this very rough 
analysis, it is impossible to make a distinction. 

Let us present again the plots of log(p) versus log(a). They are shown in Figs. 10 and 
11. We see that the case I is practically indistinguishable from a A-term, while the case II 
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is more similar to the situation in Sec. II, with different scaling regimes. It is interesting to 
note that the regimes are now three, even if the last one seems to be important only in the 
remote future. We stress that also in this case, as well as in Sec. II, w v is not constant. 

As a final result, we plot in Figs. 12 and 13 the equation of state for the scalar field. 
Now it is impossible to show an exact analytical expression, so that we only give the plots 
(in arbitrary units) in the two examined cases. It is interesting to note that in the case I, 
although so similar to the pure cosmological constant case, we nonetheless obtain a nontrivial 
plot for the equation of state. But, clearly, this point deserves further investigation. 



IV. CONCLUSIONS 

We have discussed two particular kinds of potentials which have allowed the general exact 
integration of Friedmann equations in presence of dust (ordinary and cold dark matter) and 
scalar field. This has been achieved by performing suitable transformations of variables. 
Such transformations have not been guessed by chance, but are the results of a well known 
procedure, the N other Symmetry approach [pf^,[P|pET|,pE2"| , based on an action principle. This 



was not mentioned before because it was unnecessary to the main goal of the discussion. 
Nevertheless, it now seems appropriate to stress the power of such a procedure, which allows 
to solve cosmological equations, often giving also informations on the potential and/or the 
possible coupling between the scalar field and the curvature of spacetime, without any 
limitation on the validity of the solution itself. (For details, see the literature quoted above, 
and the references therein.) 

We have seen that, with a suitable choice of integration constants for both potentials, 
it is possible to reproduce the main recent results from supernovae (initially interpreted in 
a A-model), with considerable precision, especially in the case II of the second potential. 
(The case I does not add much to what is already known in a constant A-term model.) 
What is interesting, in our opinion, is that such kinds of models can bring to a different 
evaluation of an important quantity like fi m0 . This, in a certain sense, sheds new light on 
the exponential potential of the first type, which is usually not considered as completely 
adequate for quintessence, for instance. But we have seen that, without considering a priori 
w v as a constant and having a general exact solution, something else can be learned. As a 
matter of fact, not having an almost constant makes it impossible to treat the tracker 
condition in the usual way. Also, we get appreciable values of w v , i.e., surely less than —0.7. 

Of course, all our discussion is in part still qualitative, in that we should need to make a 
more punctual analysis of observational data, and verify the best fit with the various models, 
in order to see whether and when real differences arise. Anyway, our analysis already seems 
to confirm some of the considerations made in . 



Another important point to note is that, to be realistic and cover the whole (or, at least, 
a substantially wider) range of the life of the universe, radiation (and hot dark matter) must 
be added into the game. This could allow to study the CMBR spectrum and the formation 
of structures. But it presumably destroys the possibility to integrate the system of the 
cosmological equations, leading to the necessity of using the results we established here only 
as a guide for a more complete analysis. 
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Figure Captions 



Fig. 1 - Comparison of the state equation in Eq. (36) with an approximate quadratic 
equation and with an approximate straight line. 

Fig. 2 - The straight line represents the dependence of p m versus a (p m = Da' 3 ). The 
thick curve is for the scalar field. The two tangents show the asymptotic scaling behaviour in 
the period of dominating matter and scalar field, respectively. The bullets indicate present 
time, according to a = 0.82. Units are arbitrary. 

Fig. 3 - The derivative of log(p) with respect to log(a), shows a quick transition from 
one scaling regime to another one. Units are arbitrary. 

Fig. 4 - The plot of versus time shows that is far from being constant. 

Fig. 5 - versus Q m0 . 

Fig. 6 - This plot shows the dependence of from the redshift. 

Fig. 7 - Comparison of the distance modulus 8 (derived from Eqs. (32), (33) and (34), 
and assuming fl m0 = 0.3) (continuous line) with S (derived from a A-term model, assuming 
a Q m0 = 0.37) (dots). 

Fig. 8 - Comparison of the distance modulus 5} (derived from Eqs. (68), (69) and (71) 
in case I, and assuming Q m0 = 0.3) (continuous line) with 5 (derived from a A-term model, 
assuming f2 m0 = 0.3)) [dots]. 

Fig. 9 - Comparison of the distance modulus <5(derived from Eqs. (68), (69) and (71) in 
case II, and assuming fl m0 = 0.3) (continuous line) with a (derived from a A-term model, 
assuming fl m0 = 0.3)) [dots]. 

Fig. 10 - Case I. The straight line with slope represents the dependance of log p m versus 
log a (p m = Da~ 3 ). The horizontal thick line is for the scalar field. It emulates a constant. 
The bullet indicates present time, according to r = 1.2. Units are arbitrary. 

Fig. 11 - Case II. The straight line with slope represents the dependance of logp m versus 
log a (p m = Da' 3 ). The thick line is for the scalar field. There are three approximate scaling 
regimes. The bullet indicates present time, according to r == 0.44. Units are arbitrary. 

Fig. 12 - Case I. State equation for the scalar field. Although it perfectly emulates a 
cosmological constant, the equation is nontrivial. Units are arbitrary, but coherent with 
Fig. 13. 

Fig. 13 - Case II . State equation for the scalar field. Units are arbitrary, but coherent 
with Fig. 12. 
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Fig. 1 - Comparison of the state equation in Eq. (36) with an approximate quadratic equation and with an approximate straight line. 



Fig. 2 
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log (a) 



Fig. 2 - The straight line represents the dependence of p m versus a (p m = D a~ 3 ) . The thick curve is for the scalar field. The two 
tangents show the asymptotic scaling behaviour in the period of dominating matter and scalar field, respectively. The bullets indicate 
present time, according to Tq = 0.82. Units are arbitrary. 
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Fig. 3 - The derivative of log(p) w.r.t. log(a), shows a quick transition from one scaling regime to another one. Units are arbitrary. 

Fig. 4 
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Fig. 4 - The plot of w v versus time shows that w v is far from being constant. 



Fig. 5 




Fig. 5 - WyQ versus fl mQ . 

Fig. 6 
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Fig. 6 - This plot shows the dependence of from the redshift. 
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Fig. 7 




Fig. 8 - Comparison of the distance modulus 6 (derived from Eqs. (68), (69) and (71) in case I, and assuming H m0 = 0.3) [continuous 
line] with 6 (derived from a A-term model, assuming H m0 = 0.3) [dots]. 
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Fig. 9 
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Fig. 9 - Comparison of the distance modulus 5 (derived from Eqs. (68), (69) and (71) in case II, and assuming H m o = 0.3) [continuous 
line] with S (derived from a A-term model, assuming O m Q = 0.3) [dots]. 

Fig. 10 
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Fig. 10 - Case I. The straight line with slope represents the dependance of p m versus a ( p m = Da ) . The horizontal thick line is for 
the scalar field. It emulates a constant. The bullet indicates present time, according to To = 1.2. Units are arbitrary. 



Fig. 11 
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Fig. 1 1 - Case II. The straight line with slope represents the dependance of p versus a ( p = Da ) . The thick line is for the scalar 
field. There are three approximate scaling regimes. The bullet indicates present time, according to To = 0.44. Units are arbitrary. 

Fig. 12 




Fig . 12 - Case I. State equation for the scalar field. Although it perfectly emulates a cosmological constant, 
the equation is nontrivial. Units are arbitrary, but coherent with Fig. 1 3 . 



Fig. 13 
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Fig .13 - Case II . State equation for the scalar field. Units are arbitrary, but coherent with Fig. 12. 



